Effective Dimension Transition of the Dynamics of Granular 

Materials in the Pipe. 

Akinori Awazu* 
Department of Mathematical Sciences 
Osaka Prefecture University, Sakai 599-8531 Japan. 

o 
o 
o 

<N 
•*-> 

o 

Abstract 

Phase transition of strongly excited granular materials in 2D pipe is investigated 
' numerically. By changing the ratio between width of the pipe and the height of the 

o ' 

granular bed, we observe the transition between the 1 dimensional like state and the 
^ ' actual 2 dimensional state. Moreover, it' is found I hat ihv character of the transition 

changes when the magnitude of dissipation passes through a critical value. 
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Granular materials exhibit some complex phenomena™! These complex phenomena 
strongly depend on not only the character of external forces but also that of boundary 
condition. When granular materials in slim pipes with fluid (air, water, and so on) flow 
down or are fluidized by inducing fluid from the bottom, slugging is observed@'0 . In most 
of the investigations of such slugging in a slim pipe, the system's effective dimension is 
regarded as JHll. On the other hand, two or more dimensional phenomena, bubbling and 
channeling, occur when granular materials are in a wider vessel^. The similar tendency is 
also observed for granular materials in a vibrating vesseliHl. When the width of the vessel 
is small, most of particles construct solid structure in bulk and only particles near side 
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walls flow down to the bottorma. However, in a wide vessel, most of particles are strongly 
fluidized and construct multi-roll structure like Bernard convection! From these facts, the 
behaviors of granular materials are extremely different between that in a slim vessel and 
that in a wide vessel. Then, an important problems appear; what quantity determines 'slim' 
(one dimensional system) or 'wide' (two or three dimensional system) for a given system? 
The purpose of this paper is to make a clear view on this problem. 

Here, we simulate following simplified situation!^. The system consists of two dimensional 
particles of mass 1 and diameter d in a two-dimensional box under uniform gravity. The 
width of the box is a, and the height of the box is infinite. We employ the following particle 
model which is one of the simplest model of granular materials. The equation of the motion 
of the zth particle is 

N 

X; = - 6 ( d ~ l x * - X J IK Vf/ ( d - l x * - x jD + - Vj)} + g (1) 

U(d- Ixi-Xjl) = ~(d- Ix.-x,!) 2 (2) 

here, 9 is Heviside function, N is the total number of particles, k and r\ are respectively the 
elastic constant and the viscosity coefficient^, and Xj(xj, ?/j), Vi(v x .,v V i), and g = (0, — g) 
are, respectively, the position, the velocity and the gravity of ith particles. In this model, 
the effect of particles' rotation is neglected. The system is driven by a simple energy source 
at the bottom of the box; a particle hitting the bottom with velocity (v x ,v y ) bounces back 
with the velocity (v x , V(V > 0)). We regard the bottom of box as x axis (y — 0), and x = 
as the center of box. The side walls of the box are put along x = a/2 and x = —a/2, and 
the viscosity which works between these walls and particles is zero. At the initial condition, 
we put particles bed with height b on the bottom of box. We simulate this system with 
some combinations of parameters (r),a,b), where a and b are enough large compared to the 
particle's diameter d. The above equations are calculated with the Euler's scheme. The 
time step St is set enough small such that 5x, the displacement of the zth particle during St, 
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does not exceed a given value. In this paper, we set (—g,V) so that the average height of 
the center of mass of the system CM =< (J2Ui)/N > t keeps enough large compared to b. 

Figure 1 shows typical snapshots of the system for respectively, (a) a « b with 77 > 77*, 
(b) a >> b with 77 > 77*, (c) a << b with 77 < 77*, and (d) a >> b with 77 < 77*. Here, we fixed 
6, whereas (— g, V) of (a) is same as that of (b), and (— g, V) of (c) is same as that of (d). In 
cases of (a) and (c), only the particles distribution in the horizontal direction is symmetric, 
and the center of mass of this system moves a little only in the vertical direction. On the 
other hand, the particles distribution is non-uniform in vertical and horizontal direction, 
and one convection appears in cases of (b) and (d). In order to characterize the system, we 
introduce the order parameter L(t) = J2i \xi{t)v yi (t)\/N which indicates the strength of the 
convection of the system. Figure 2 (a) and (b) are typical probability distributions of L(t) 
which is given by one time series for respectively the cases of a << b and those of a » b. In 
Fig. 2 (a), the peak of the probability distribution of L(t) appears at L(t) = which means 
there are no convection for the case a « b. In these cases, the effective dimension of this 
system can be regarded as 1. We name such states as the 'ID state'. On the contrary the 
peak of the probability distribution of L(t) appears at L(t) > in Fig. 2 (b) which means 
a convection with finite magnitude appears in the system. In these cases, the system is 
actually 2 dimensional system. We name such states with a convention as the '2D state'. 
Such a transition between the ID state and the 2D state which depends on the relation 
between a and b is observed in a simple system. 

In our simulation, particles are strongly excited by the energy source at the bottom of the 
box. Then, if we set a improper g for a small 77 or a small b, the height of particles diverge. 
In order to clear off such cases, we need to control g for several 77 and b. By the simulation, 
we found the fact that CM have almost same values independent of a with a « b for a 
fixed set of g, 77, b and V. Hereafter, we fix V and determine the gravity g for a given set 
(77, b) independent of a. In this paper, we set V = 2.0 and g(rj, b) with which CM ~ 6O0? is 
realized for all cases of a « b. When we set V = 3.0 or V = 4.0, we can get qualitatively 
same results as those of following discussions with V = 2.0. 
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Before the discussion of the transition width at which the 1D-2D transition takes place 
or the character of the transition, we discuss the dependency of 77 for the ID states. Figure 
3 (a) shows the time averaged packing fraction profile of the y direction as the function of 

V — Vcm for several r\ with a = 9d and b = 30d. Here, ycM is the y component of center of 
masses of particles, and we define the packing fraction as followings. We divide the space by 
a lattice with the lattice constant d, and we define the number of centers of particle within 
each d x d square as the packing fraction in the square. The packing fraction of y direction 
is given as the average of them through the x direction for each y. For large 77 (77 = 0.6, 1.0), 
each profile includes a flat region with large packing fraction near the center of mass of 
this system. These flat regions mean the existence of the solid structure in which particles 
are almost completely packed. The length of this flat region decreases with decreasing 77, 
and this length becomes for 77 = 77* ~ 0.38. If 77 < 77* (77 = 0.25,0.35), each packing 
fraction profile includes no flat region, and maximum packing fraction is smaller than that 
of 77 > 77*. It means that no solid structures are created for such small 77. Thus, a transition 
between a state which includes a solid structure and the other state which include no solid 
structure occurs at the critical value 77 = 77*. Similar results are obtained in following two 
cases, b = 20d and b = 40d. Now we introduce following no dimensional values: a' = a/d, 

V = b/d, and e' — 1 — e where e = exp(irr]/(k — f] 2 )^)), and the length of flat regions 
h(e') = h'(e')d. (h'(e') has no dimensions.) Here, e indicates the coefficient of restitution 
for head-on collisions between two particles!. The relation between two rescaled values, e'b' 
and h y 2 ~ ; is obtained as shown in Fig. 4(a). The rescaled critical point (e'b')* is determined 
independent of the system size, and the profile of this relation fits with 

^ = h' \e'b' - (e'b'yr (3) 

near (e'b')* ~ 0.282. (h' Q is constant.) Moreover, we introduce the cluster length u(e') = 
a'(e')d (cr'(e') has no dimensions.) which is defined as the distance between two nearest 
inflection points from the maximum point of the packing fraction profile. Thus, the relation 
between two rescaled values, e'b' and a'(e')/b', are obtained as shown in Fig.4 (b). The 



profile of this relation fits with 
a'(e') 0.0255 0.105 n nn 

= w ~ ^ + °' 98 ' (4) 

We expect a'(e')/b' -> oo for e'V -> 0, and a'(e')/b' ~ 1 for e'V -> oo. 

Now, we discuss the character of transition between the ID state and the 2D state. 
Here, we define this transition width a* as the maximum width that the 2D state cannot 
be observed in the system. Figure 5 shows the typical probability distributions of L(t) for 
several width a around a = cr(e'). Here, (a) indicates for the case e' < e'* (rj < rj*) and 
(b) indicates for the case e' > e'* (rj > ?]*)both for b = 20d. In Fig. 5 (a), each profile of 
probability distribution of L(t) includes only one peak. The position of peak is at L(t) = 
for a < cr(e'). For a > cr(e'), however, the peak appears at L(t) > and this peak moves with 
the width of the system. Thus the continuous transition between the ID state and the 2D 
state appears for the case e' < e'* (rj < rf), and the transition width is given as a* ~ cr(e'). 
In Fig. 5 (b), on the contrary, probability distributions of L(t) for a, which is a little larger 
than cr(e'), include two peaks at L(t) = and L(t) > 0. This means that two locally stable 
states, one is the ID state and the other is the 2D state, coexist and they appear periodically 
for the case of e' > e'* (rj > rf). In such cases, the system includes a solid structure when 
ID state is realized. In this solid structure, the friction between particles are strong because 
particles are densely packed. Hence, the solid structure is break-proof, and this originates 
the stability of the ID state for a > cr(e'). Figure 6 shows the semi-log scale profiles of 
probability distributions of L(t) for respectively a < cr(e'), a ~ cr(e') and a > cr(e'). When 
a' < cr'(e'), the profile is proportion to exp(— r yoL 2 ), which means that fluctuations of L{t) 
are so small that they can be neglected. However, a profile which proportion to exp(— jiL) 
is obtained when a is close to cr(e'), which means that the fluctuation from L(t) = become 
large. Moreover, the profile includes two peaks at L(t) = and L(t) > when a > cr(e'). 
Then, the transition width a* is regarded as a* ~ cr(e') for also e' > e'* (rj > rf). These 
results mean following two facts. I) The critical width which the 2D state can appear, 
is equal to the cluster length. II) The magnitude of dissipation separates the type of the 



transition between the ID state or the 2D state. We obtain the similar results for the case 
V = 30. Then, by using rescaled parameter e'b' and a'/b', the phase diagram is obtained as 
shown in Fig. 7. 

In this paper, we simulated strongly excited granular materials in a 2 dimensional pipe. 
When the width of the system is small, a characteristic state, in which the effective di- 
mension of particles' dynamics is regarded as 1, appears. We named this state as the ID 
state. However, when the width of the system is larger than a critical width, the actual 2 
dimensional state in which convection appears. We named this state as the 2D state. More- 
over, we found that the critical width is equal to the length of the cluster which appears 
in the ID state. The character of the ID state depends on the magnitude of dissipation as 
followings. When the magnitude of dissipation is larger than the critical value, the system 
include solid structure. On the contrary, such solid structure disappears when the magni- 
tude of dissipation is smaller than the critical value. According to the differences of the 
magnitude of dissipation, following two types of behaviors appears in the system near the 
critical width. When the magnitude of dissipation is smaller than the critical value, the 
continuous transition between the ID state and the 2D state appears. On the contrary, 
when the magnitude of dissipation is larger than the critical value, two meta-stable states, 
the ID state and the 2D state, appear periodically for a little over the critical width. By 
another simulation, the existence of such a critical magnitude of dissipation is reportedSE!. 
The analytical derivation of this critical value is one of the most important issue for the 
research of granular materials. Moreover, simulations of more highly excited systems, larger 
systems, and analytical study of the critical width for several magnitude of dissipation are 
important future issues. 

The author is grateful to H.Nishimori for useful discussions. This research was supported 
in part by the Ibaraki University SVBL and Grant-in-Aid for JSPS Felows 10376. 
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FIGURES 

Fig. 1. Illustration of 2D pipe and typical snapshots respectively, (a) a « 6 with rj > 77*, (b) 
a >> b with 77 > 77*, (c) a « b with 77 < 77*, and (d) a >> b with 7/ < 77*. 

Fig. 2. Probability distributions of L(£) respectively (a) a << 6, and (b) a >> 6. 

Fig. 3. Packing fraction profile for y direction respectively i] = 0.25,0.35,0.6, and 1.0. 

Fig. 4. Relation between two rescaled values, (a) e'b' and h ^ - , and (b) e'6' and a'(e')/b' 

Fig. 5. Typical probability distributions of L(t) of (a) e' < e'* (77 < 77*) and (b) e' > e'* (77 > 77*) 
for several width around a critical value a = cr(e'). 

Fig. 6. Semi-log scale profiles of probability distributions of L(t) for respectively a < cr(e'), a is 
a little smaller than er(e') and a > cr(e') with e' ~ 0.247 (7/ = 1.0) 

Fig. 7. Phase diagram of the effective dimension of the system which depends on e'6' and a'/b' 
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